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This  paper  which  is  based  on  neighboring  optimal  control  investigates  the  application  of  an 
analytical  feedback  guidance  scheme  for  the  problems  of  low-thrust  orbital  transfer.  Neigh¬ 
boring  optimal  control  is  one  of  the  methods  that  has  been  proposed  for  the  closed-loop 
optimal  guidance.  This  method  focuses  on  expanding  the  cost  function  to  second  order  and 
examining  the  robustness  of  the  system  by  the  feedback  law  gains.  At  first,  the  open-loop 
analytical  guidance  policy  is  considered  as  the  optimal  thrust  steering  program  that  will 
transfer  the  vehicle  from  the  inclined  low  earth  orbits  to  the  high  earth  orbits.  Secondly, 
proper  feedback  optimal  guidance  laws  are  analytically  obtained  to  maintain  the  trajecto¬ 
ries  around  their  optimum  with  assessed  disturbances.  The  proposed  guidance  scheme  is 
distinguished  for  two  desired  performance  indices  as  minimum-time  and  minimum-effort. 
Finally  performance  indices  are  compared  and  the  best  policy  is  obtained  regarding  the  ro¬ 
bustness  of  the  two  performance  indices  against  disturbances. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 


Optimal  control  of  dynamic  systems  can  be  divided  into  two  groups  of  open-loop  and  closed-loop  optimal  controls. 
Open-loop  optimal  control  is  time-dependent  and  initiates  from  original  conditions,  but  closed-loop  optimal  control  is  state- 
dependent  variables  and  has  robustness  against  disturbances.  In  this  way,  many  methods  can  be  illustrated  in  the  field  of 
applied  optimal  control  as  dynamic  programming,  calculus  of  variations,  neighboring  optimal  control,  Pontryagin’s  min  (or 
max)  principle,  etc. 

Dynamic  programming  (DP)  requires  optimal  control  over  an  n-dimensional  grid  near  the  optimal  path.  This  method 
needs  storing  a  huge  amount  of  data;  by  contrast,  another  method  as  neighboring  optimal  control  (NOC)  needs  only  storing 
nominal  path  and  control.  Actually,  neighboring  optimal  can  be  generated  as  closed-loop  optimal  control  by  feed-forward 
of  the  control  as  nominal  optimum  plus  feedback  of  the  deviation.  The  fundamental  of  NOC  is  based  on  minimizing  the 
second-order  terms  in  the  expansion  of  the  performance  index  function  about  an  optimal  nominal  trajectory  [1], 

Based  on  modification  of  a  well  known  neighboring  optimum  feedback  control  scheme  and  shooting  method,  Pesch  pro¬ 
posed  neighboring  optimum  guidance  of  a  space  shuttle  orbiter  in  Ref.  [2],  Next,  NOC  with  respect  to  a  feedback  law  for  the 
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Nomenclature 

NOC  Neighboring  optimal  control 

x(t)  States  of  the  system  equation 

D(t)  Vector  of  controls  in  general 

x(t)  Time  derivative  of  states 

J  Performance  index 

tf  Final  time 

L  Operator  of  performance  index 

H(t)  Hamiltonian  function 

A.(t)  Co-state 

|  Magnitude  of  acceleration 

P  (t)  Control  variable  for  low-thrust 

i(t)  Orbital  inclination 

V(t)  Orbital  velocity 

a(t)  Semi  major  axes 

fjt  Gravitational  constant 

Su  Variation  of  control  vector 

K(t)  NOC  gain 

Sx  Variation  of  state  vector 

T  Division  operator  for  NOC 


advanced  system  was  investigated  by  a  robust  feedback  algorithm  for  a  near  minimum-fuel,  two  stages  launch  vehicle,  by 
Seywald  et  al.  [3].  A  neighboring  optimal  adaptive  method  has  been  proposed  by  Balakrishnan.  He  used  neighboring  optimal 
control  and  linear  optimal  guidance  law  for  his  neural  network  approach  [4],  An  aero-assist  optimal  orbital  transfer  with 
fuel  consumption  was  devised  by  a  neighboring  optimal  guidance  scheme  for  a  nonlinear  dynamic  system  [5]. 

Schaub  and  Alfriend  developed  an  impulsive  feedback  control  to  establish  specific  relative  orbits  for  spacecraft  formation 
flying  [6],  In  Ref.  [7],  the  technique  of  neighboring  optimal  control  was  proposed  to  compute  near  optimal  trajectories.  In 
their  work,  Jardin  et  al.  computed  trajectories  in  general  wind  fields  with  nominal  solutions  based  on  the  Zermelo  problem, 
and  the  excellent  performance  of  this  method  was  investigated.  Also,  in  Ref.  [8]  Jardin  compared  neighboring  optimal  control 
to  LQR  for  commercial  airlines. 

Ref.  [9]  presented  near  optimal  control  laws  for  minimum-fuel  rendezvous  between  satellites  of  arbitrary  eccentricity 
in  elliptic  orbits.  Axelrod  and  Gualman  in  Ref.  [10]  designed  a  feedback  law  that  guided  the  spacecraft  along  an  optimal 
neighboring  trajectory.  Their  work  was  based  on  electrical  propulsion  with  discrete  thrust  levels.  Gros  et  al.  proposed  an 
approach  by  considering  neighboring  extermals  for  multiple  inputs  with  a  feedback  law  [11], 

Hamel  and  Lafontaine  in  Ref.  [12]  developed  a  feedback  control  law  for  neighboring  fuel  optimal  reconfiguration  of  for¬ 
mation  flying  spacecraft.  Wurth  and  Hannemann  used  neighboring  extermals  to  solve  dynamic  optimization  for  suitable 
application  in  a  nonlinear  model-predictive  control  and  dynamic  real-time  optimization  [13],  Jiang  worked  on  methods 
to  achieve  trajectory  with  neighboring  optimal  control  for  reusable  launch  vehicle  [14],  Liu  et  al.  investigated  an  optimal 
control  for  periodic  tasks  of  discontinuous  dynamic  systems.  They  used  differential  dynamic  programming  to  optimize  and 
generate  an  optimal  control  law  regarding  a  linear  local  model  in  the  neighborhood  of  the  nominal  trajectories  [15].  In 
Ref.  [16]  neighboring  external  regarding  to  use  error  compensation  is  investigated  by  Kim  and  Hall. 

The  optimal  low-thrust  orbit  transfers  have  received  a  great  amount  of  attention  in  flight  mechanics  and  astrodynamics 
literatures.  The  evolutions  of  low-thrust  propulsion  technologies  have  reached  a  point  where  such  systems  have  become 
an  economical  option  for  many  space  missions.  Also  the  developments  of  efficient  control  laws  have  received  an  increasing 
amount  of  attention  in  recent  years,  and  few  studies  have  examined  the  subject  of  inclination  changing  maneuvers  [17],  In 
Ref.  [  18]  the  original  Edelbaum  theory  for  low-thrusters  was  revisited  within  the  framework  of  optimal  control  for  min-time. 
In  Ref.  [19],  Casalino  investigated  Edelbaum’s  approach  again,  but  some  improvements  were  introduced,  while  maintaining 
the  assumption  of  quasi-circular  orbits.  Shafieenejad  and  Novinzadeh  in  2010  had  more  research  works  to  achieve  analytical 
open-loop  solutions  for  low-thrust  orbital  maneuvers.  They  examined  and  compared  different  performance  indices  for  low- 
thrust  optimal  control  problems  [20,21], 

The  current  research  augments  the  previous  research  works  [20,21]  as  studying  analytical  closed-loop  guidance  regard¬ 
ing  neighboring  optimal  control.  The  resulting  optimal  guidance  law  is  composed  of  the  analytical  open-loop  controls  (based 
on  Shafieenejad’s  research  works  [20,21])  and  analytical  closed-loop  time  varying  guidance  law  with  neighboring  optimal 
control  (NOC).  Hence,  in  the  middle  section  of  this  work,  the  open-loop  optimal  low-thrust  orbit  transfers  will  be  intro¬ 
duced  with  respect  to  the  previous  works  and  the  NOC  method  will  be  described.  Then,  the  NOC  solutions  to  the  low-thrust 
as  minimum-time  and  minimum-effort  will  be  investigated.  This  work  concludes  many  different  simulated  scenarios  to 
show  the  robustness  of  the  NOC  for  two  performance  indices.  Finally,  comparisons  between  these  two  performance  indices 
will  be  made  with  respect  to  their  ability  to  damp  disturbances  that  are  exerted  on  the  system. 
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Fig.  1.  Low-thrust  orbital  transfer. 


2.  Open-loop  optimal  control  for  low-thrust  orbital  transfer 

The  calculus  of  variation  is  utilized  in  this  paper  and  formulated  briefly.  At  first  one  can  consider  the  mathematical 
representations  of  this  system  as 


x(t)  =  f[x(t),  u(t),  t] ,  t0  <  t  <  tf 


(1) 


u(t)  and  x(t)  are  considered  as  the  vector  of  m-control  component  and  n-state  of  the  system.  The  conventional  form  of  the 
cost  function  is  expressed  by  Eq.  (2)  without  considering  the  terminal  part: 


(2) 


In  the  calculus  of  variations,  the  Hamiltonian  of  the  system  is  determined  as  follows: 


H  =  L(x,  u,  t)  +  XT(t)f(x,  u,  t). 


(3) 


A.(t)  demonstrate  the  co-state  of  the  system.  In  the  next  step,  the  optimality  condition  of  the  system  is  considered  as  the 
main  principle  of  optimality  (3 H/du)  =  0.  All  mentioned  equations  need  to  be  satisfied  simultaneously  by  the  boundary 
conditions  of  the  optimal  control  problem.  From  optimality  conditions,  m-control  components  are  achieved  as  functions  of 
x(t)  and  X(t).  In  the  most  practical  applications  in  optimal  control,  we  confront  difficulties  of  nonlinear  equations.  In  this 
way,  it  is  very  important  to  facilitate  the  problem  and  overcome  difficulties  of  analytical  solutions  [1,19,20]. 

3.  Low-thrust  orbital  transfers 

In  this  work,  problems  of  low-thrust  orbital  transfers  are  formulated  in  the  framework  of  optimal  control  problems.  In 
this  formulation  the  thrust  direction  is  considered  as  an  optimal  control  variable.  The  advantages  of  these  guidance  laws 
are  that  they  are  the  exact  solutions  of  the  two-point  boundary  value  problems  that  can  be  used  later  for  the  analytical 
closed-loop  optimal  guidance  laws.  Also,  terminal  conditions  of  time-free  plane  changing  maneuvers  of  low-thrust  space¬ 
crafts  are  satisfied.  Furthermore,  in  this  research,  the  several  difficulties  that  are  associated  with  the  numerical  determi¬ 
nation  of  optimal  open-loop  and  closed-loop  solutions  for  the  nonlinear  systems  are  eliminated  by  proposing  analytical 
solutions. 

Edelbaum  developed  simple  and  applicable  equations  by  using  planetary  equations.  In  this  way,  an  analytical  optimal 
control  solution  is  important  for  space  researchers.  It  should  be  noted  that  analytical  solutions  satisfy  boundary  conditions 
and  performance  measures.  There  are  some  suppositions  in  the  Edelbaum  equations.  The  pitch  angle  is  considered  as 
a  control  variable;  orbital  inclination  and  orbital  velocity  are  also  considered  as  state  variables.  The  acceleration  of  the 
spacecraft  is  constant  and  each  orbital  revolution  is  a  near-circular  orbit  regarding  the  Edelbaum  analysis  [  17], 

Minimum-time  and  minimum-effort  performance  indices  are  studied  again  for  non-mass  variant,  and  the  results  are 
examined  for  many  maneuvers  from  LEO  to  GEO.  From  the  schematic  diagram  of  Fig.  1,  f}(t)  is  the  out-of-plane  or  thrust 
angle  and  f  represents  the  magnitude  of  the  acceleration  vector.  Two  components  £t,  th  are  considered  as  =  I  cos  (fi) 
and  fh  =  ?  sin  (fi).  |t,  f/,  represent  the  tangential  and  normal  components  of  the  acceleration  vector  respectively.  Equations 


I.  Shafieenejad  et  al.  /  Mathematical  and  Computer  Modelling  58  (2013)  1714-1726 


1717 


Fig.  2.  Optimal  orbital  maneuver  from  LEO  to  GEO. 


of  this  orbital  maneuver  for  low-thrusters  are  as  follows: 
di  _  2|sin0 0) 
dt  ~  ttV 
dV 

|-=-|cos(«. 

The  relationship  between  orbital  velocity  and  semi  major  axes  can  be  achieved  by  the  energy  equation  ^ 

[17-20], 

3.1.  Min-time  low-thrust  orbital  transfers 

In  this  part,  the  transfer  problem  is  cast  as  a  min-time  problem  between  initial  and  final  conditions  i0,  V0  and  if,  Ir¬ 
respectively.  Co-states  A.,,  A,„  are  achieved  analytically  as  functions  of  optimal  control  respectively:  A.,-  =  ^  7tVsf'^\lv  = 
(see  Refs.  [20,21]). 

The  time  derivatives  appearing  in  the  governing  equation  can  be  written  with  respect  to  fi.  In  this  way,  now  fi  becomes 
the  implicit  independent  variable.  ^  and  the  new  system  equations  are  given  by 

di  _  2 

(5) 

dV  _  V  cos  (/l)  1  ; 

dp~  sinOS)  ' 

Due  to  the  system  Eq.  (5),  the  results  as  a  function  of  the  control  angle  P  are  achieved: 
sinCfr)^ 
sin  (P) 

,  =  ll  +  2AL-M 

Obviously,  it  is  necessary  to  specify  the  values  of  p0  and  fif.  This  can  be  accomplished  through  using  the  known  initial  and 
terminal  conditions  and  solving  a  set  of  nonlinear  algebraic  equations.  As  mentioned,  solutions  have  been  expressed  in  terms 
of  P  and  it  is  important  to  get  an  analytical  closed-loop  guidance  law  based  on  neighboring  optimal  control.  Fig.  2  shows 
how  optimal  orbital  maneuver  is  accomplished  around  the  Earth.  For  the  next  step,  time-to-go  is  derived  with  respect  to  p 


(6) 

(7) 
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as  follows: 


Time  to  go 


VfSm(Pf-P) 
I  sinG0) 


(8) 


3.2.  Min-effort  low-thrust  orbital  transfers 


In  this  part,  after  considering  minimum-time,  the  transfer  problem  is  cast  as  minimum-effort  and  as  in  the  last  section 
i,  V  are  considered  as  the  state  variables.  The  performance  index  is  considered  as  J  =  p  (t)2dt.  Also  in  this  part,  co-states 
are  given  for  this  performance  index  as  A,-  =  ^  n  ;  ZcosW> ,  ^  _  p(cas(p)p-2sm(p)) .  see  Refs  [20,21]. 

As  mentioned  in  the  previous  section,  the  time  derivatives  can  be  derived  again  with  respect  to  p.  On  the  other  hand  an 
implicit  relation  of  control  angle  and  time  is  obtained;  therefore  the  new  system  equation  is 

di  2  sin(/0)  (2  +  /S2) 

W  ~  ;ry8(sinG0)/0  +  2cosG0)) 
dV  V  cos 06)  (2  +  p2) 

dp  ~~p(sin(p)p +  2  cos  (p))' 

If  time  is  considered  as  an  independent  variable,  exact  closed-form  solutions  that  are  suitable  for  neighboring  optimal  control 
as  analytical  feedback  guidance  law  are  inaccessible.  Hence,  state  variables  of  the  system  equations  as  functions  of  control 
angle  help  us  in  the  next  parts  to  achieve  NOC: 

VfPf{sin(Pf)Pf  +  2cos(Pf)) 

P  (sin(P)P  +  2  cos(P)) 

P  (sin(P)P  +  2  cos  (P))  jo 

- VJtTW) - 

Through  using  the  known  initial  and  final  conditions  and  solving  a  set  of  nonlinear  algebraic  equations 
obtained.  Time-to-go,  according  to  p,  is  derived  for  this  performance  index: 


(10) 

(11) 

,  p0  and  Pf  can  be 


Time  to  go 


rf>f  /  Vf  Pf  (2  +  p2)  (sin(P)p  +  2  cos  (P))  \ 

Ip  \P2I  (-4/8  sin (P)  cos  (P)  -  4  cos  (P)2  -  P2  +  P2  cos  (P)2)  J 


(12) 


4.  Closed-loop  optimal  control  via  the  neighboring  optimal  control  (NOC)  method 

One  of  the  main  problems  that  plays  an  important  role  in  space-systems’  motion  is  determining  the  optimal  trajectory  and 
providing  a  proper  closed-loop  optimal  guidance  law  to  maintain  the  system  while  disturbances  are  being  assessed.  For  the 
first  step  to  achieve  this  purpose,  the  solution  of  the  desired  problem  based  on  the  open-loop  optimal  control  theory  should 
be  determined.  In  this  way,  the  calculus  of  variation  is  one  of  the  most  important  methods  to  provide  proper  solutions  for 
optimal  control  problems.  This  method  takes  even  more  credit  when  the  problem  has  a  nonlinear  behavior  as  the  presented 
problem  in  this  research.  Providing  an  analytical  optimal  solution  is  very  difficult  due  to  complexities  of  equations  and 
advanced  mathematical  methods;  however,  achieving  an  analytical  optimal  solution  for  the  desired  problem  is  very  useful 
and  it  can  assist  us  to  achieve  analytical  closed-loop  and  enhance  the  robustness  of  the  system  against  disturbances. 

The  purpose  of  this  part  is  providing  a  modified  optimal  control  vector  to  damp  the  exerted  disturbances  based  on  an 
analytical  neighboring  closed-loop  control  method.  This  method  which  results  in  the  time-variant  feedback  control  law 
minimizes  the  second-order  optimality  criterion  of  the  nominal  trajectory  disturbances.  According  to  the  NOC  method,  the 
gain  feedback  can  be  planed  as  a  time-interval  function  between  the  end  time  and  the  current  time,  which  predicts  the 
variation  of  final  time.  One  of  the  advantages  of  this  method,  unlike  other  methods  such  as  LQR,  is  the  possibility  to  modify 
the  end  time.  It  should  be  noted  that  the  NOC  method  is  based  on  a  feedback  algorithm  which  means  that  the  control  variable 
can  be  reproduced  and  rectified  by  Su  [1,7]: 

<$u(t)  =  Ku(t)<5x(t)  (13) 

where  <5u(t)  represents  the  variation  of  the  control  vector  which  generated  by  the  gains  and  <5x  denotes  the  disturbed  state 
vector.  Also  I<(t)  is  the  neighboring  time-variant  optimal  feedback  gain. 

Since  the  end  time  of  the  problem  in  calculations  is  not  specified,  the  technique  presented  in  this  paper  can  be  employed 
to  calculate  the  end-time  variations  based  on  the  gains  and  rectify  the  nominal  optimal  response  subsequently: 

dtf  =  K  t(t)Sx(t). 


(14) 
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Fig.  3.  Schematic  plan  of  the  closed-loop  law  based  on  the  NOC. 


According  to  this  method,  the  feedback  rule  for  p  and  f);  as  control  variable  and  the  value  of  the  control  at  end  time 
respectively  are  expressed  as  the  below  matrix  equation: 


(«) 


(£)■ 


(15) 


So  the  disturbances  exerted  to  variables  such  as  X  are  defined  as  follows: 


*(*)  =  (*-*-)■ 

On  the  other  hand,  the  partial  differential  of  P  and  Pf  to  state  variables  are 


3  P  , 

(1\ 

9x2 

9x,  1 

l  r) 

JPf 

dp  , 

(~V 

\  9xi 

9x2  1 

[r. 

)Wf 

dPf 

(~ 1 

\  9x2 

dxt 

v  r 

)Jp 

dpL  = 

(-) 

9xi 

9x2 

3 T 

The  presented  division  operator  T  in  the  above  equations  is  also  defined  as  follows: 


(16) 

(17) 

(18) 

(19) 

(20) 


'-(SXSMSKB)- 

Finally,  control  variables  with  regard  to  disturbances  are  defined  as  (see  Fig.  3) 

Unew  =  U0pt  +  SU.  (22) 

And  in  the  next  step,  final  time  variation  according  to  the  following  equation,  where  T  is  time  as  a  function  of  P,  fy,  is 
demonstrated  [1,7]: 


dtf 


17  £ 

9xi 


*L*k\,x,(VL*!L 

dpf  dxj  \  dp  9x2 


d^dpA 
dpf  9x2  ) 


(23) 


Calculating  the  above  equation  necessitates  achieving  final  time  and  states  with  respect  to  P,  Pf. 

In  the  present  paper  an  analytical  method  is  employed  to  get  feedback  gain;  however,  other  methods  like  backward 
sweep  integration  that  is  presented  in  Ref.  [1]  confront  researchers  with  computational  problems.  But  due  to  analytical 
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Fig.  4.  NOC  response  of  the  orbital  inclination  based  on  disturbance  of  type  (I)  for  min-time. 


Fig.  5.  NOC  response  of  the  orbital  velocity  based  on  disturbance  of  type  (I)  for  min-time. 


solutions  provided  for  this  problem,  differentiation  regarding  no-singularity-existence  along  the  time  domain  is  easier  and 
more  suitable.  It  should  be  noted  that  this  algorithm  is  applicable  for  problems  that  the  end  time  is  not  fixed. 

According  to  the  method  proposed  by  the  optimal  neighboring  method,  the  desired  gains  will  be  achieved  and  we  can 
control  the  space  system  about  the  near  optimal  trajectory  by  multiplying  those  gains  by  the  deviation  of  the  state  variables: 

Su  =  Ki  (t)  Sxi  +  K2  (t)  8x2  .  (24) 

As  shown  in  Fig.  3,  the  new  state  variables  are  obtained  through  applying  disturbances  to  a  system  and  measuring  the  re¬ 
sulted  values.  Subsequently  the  new  gain  which  leads  us  to  compute  the  offset  values  of  the  control  variable  or  8u  is  achieved 
by  comparing  the  new  values  with  the  optimum  values  of  state  variables.  The  summation  of  this  value  and  the  basic  uopt 
value  which  is  obtained  through  the  open-loop  optimal  control,  and  considering  the  deviation  of  state  variables  from  the 
open-loop  optimal,  results  in  a  new  control  variable  unew. 

4.  f .  Closed-loop  optimal  control  of  a  low-thrust  spacecraft  orbit  transfer  with  min-time  criterion 

As  described  before,  open-loop  optimal  control  results  based  on  min-time  criterion  as  a  function  of  fi  and  fy  are 
considered  primarily  ( ft  is  the  control  variable ).  Considering  Eq.  (15),  <5/3  and  <5/3/  will  be  achieved  by  applying  the  T  operator 
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Fig.  7.  NOC  response  of  the  orbital  inclination  based  on  disturbance  of  type  (II)  for  min-time. 


as  defined  in  Eqs.  (17)-(21).  In  the  next  step,  the  value  of  dtf  based  on  Eq.  (23)  is  achieved  and  the  new  control  variable 
5>S(t)  =  /Snomi  +  is  resulted  according  to  the  defined  algorithm  (Fig.  3). 

The  equations  that  are  presented  in  this  paper  for  analytical  closed-loop  clearly  indicate  that  the  end  time  is  not  fixed. 
Also,  it  is  necessary  to  rectify  the  values  of  final  time  ts  and  the  final  angle  fi;  that  are  presented  in  Table  1.  It  should  be  noted 
that  calculating  new  y3/,  tf  is  one  of  the  advantages  of  the  NOC  method  against  other  methods  to  get  feedback  gains: 


a0  =  7000  (km)  a/  =  42,000  (km)  i0  =  28.5° 


if  =  0°  $  =  3.5  x  10 


With  respect  to  Ref.  [16],  to  verify  the  closed-loop  performance  based  on  the  neighboring  method,  it  is  proper  to  apply 
deviations  and  disturbances  to  the  system  and  try  to  guide  the  system  about  the  near  optimal  neighborhood.  Therefore, 
in  this  part,  two  types  of  disturbances,  as  initial  deviations  and  noises,  are  applied  to  the  system  in  order  to  study  the  ro¬ 
bustness  of  the  system.  Disturbances  of  type  (I)  are  modeled  as  initial  deviation  with  amplitude  of  0.01  values  of  states  as 
Si0  =  [0.01]  i0,  8V0  =  [0.01]  V0-  In  type  (II)  a  sinusoidal  disturbance  which  has  an  amplitude  of  0.05  of  control  variable 
is  exerted  on  the  control  variable  to  investigate  the  effect  of  NOC  gains.  State  variables  are  achieved  in  the  next  figures  for 
three  different  manners:  the  open-loop  optimal  solution  without  disturbance,  the  open-loop  optimal  with  disturbance  and 
neighboring  optimal  control  solutions  (NOC). 
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Table  2 

Open-loop  and  closed-loop  optimal  control  results  based  on  the  NOC  method  for  min-effort. 


Solutions 

Min-control-effort 

Open-loop  optimal  control 

Open-loop  final  time  (day) 

191.3755 

PoO 

23.8226 

fir  O 

61.5515 

Ap  control  command  (°) 

37.72887 

Closed-loop  feedback  guidance  law  (NOC) 

NOC  final  time  (day)  (disturbance  type  I) 

191.3769 

NOC  Pf  (°)  (disturbance  type  I) 

61.6194 

NOC  final  time  (day)  (disturbance  type  II) 

191.4078 

NOC  Pf  (°)  (disturbance  type  II) 

61.6236 

Fig.  11.  NOC  response  of  the  orbital  velocity  based  on  disturbance  of  type  (I)  for  min-effort. 

As  shown  in  Figs.  4-9,  the  presented  guidance  law  can  guide  the  spacecraft  about  the  neighbored  of  the  open-loop  optimal 
trajectory. 
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Fig.  13.  NOC  response  of  the  orbital  inclination  based  on  disturbance  of  type  (II)  for  min-effort. 


4.2.  Closed-loop  optimal  control  of  low-thrust  orbital  transfer  with  the  min-effort  criterion 


In  this  section,  the  analytical  closed-loop  of  the  orbital  transfer  problem  based  on  a  min-effort  criterion  is  investigated. 
As  the  first  step  according  to  the  open-loop  optimal  solutions  of  the  previous  part  (see  Eqs.  (10),  (11)),  the  state  variables 
as  a  function  of  fi  and  fis  are  rewritten.  Then  the  feedback  gains  are  achieved  by  considering  the  NOC  closed-loop  Eqs. 
(15)-(21).  As  shown  in  Fig.  3,  disturbance  is  exerting  on  the  system  and  its  effect  on  the  state  variable  is  considered.  Also 
Figs.  10-15  demonstrate  that  the  open-loop  solution  has  high  sensitivity  to  disturbance  and  a  satellite  tends  to  deviate 
from  the  optimal  trajectory,  but  in  the  closed-loop  control  solution,  all  disturbances  have  been  damped  perfectly.  Hence, 
the  optimal  trajectory  is  very  close  to  the  optimal  solution  without  disturbances.  To  study  the  effect  of  disturbance  on  the 
state  variables,  the  following  disturbances  in  two  types  are  added  to  the  system.  Disturbance  of  type  (I)  is  modeled  as  a 
sinusoidal  disturbance  with  amplitude  value  of  0.1.  In  type  (II),  a  random  disturbance  which  has  a  normal  distribution,  unit 
variance  and  zero  mean  value  is  exerted.  Finally,  boundary  conditions  are  the  same  as  what  are  stated  in  Table  2  and  also, 
the  variation  of  fi  and  fy  regarding  the  applied  disturbances  to  the  state  variables  are  taken  into  consideration: 


a0  =  7000  (km)  af  =  42,000  (km)  i0  =  28.5° 


if  =  0°  |  =  3.5  x  10" 


■’(f)- 
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Fig.  14.  NOC  response  of  the  orbital  velocity  based  on  disturbance  of  type  (II)  for  min-effort. 


5.  Conclusions 

In  this  research,  optimal  control  problems  for  two  performance  indices  as  min-time  and  min-effort  are  introduced. 
Secondly  these  two  indices  for  low-thrust  orbital  maneuvers  are  investigated  to  achieve  analytical  closed-loop  guidance 
laws  by  Neighboring  Optimal  Control  (NOC).  This  method  needs  only  storing  nominal  path  and  optimal  control.  By  feed¬ 
forward  of  the  control  as  nominal  optimum  plus  feedback  of  the  deviation,  neighboring  optimal  can  be  generated  as  closed- 
loop  optimal  control.  It  is  concluded  from  Tables  1,  2  that  the  min-effort  criterion  is  more  appropriate  than  the  min-time 
criterion  regarding  A  f)  and  final  time.  The  results  of  this  work  demonstrate  the  fact  that  the  closed-loop  optimal  response 
has  a  better  effect  on  the  min-effort  and  it  can  damp  the  disturbance  with  higher  amplitude  for  this  criterion.  In  other  words, 
it  is  one  of  the  novelties  of  this  work  to  assign  a  better  performance  index  than  min-time  regarding  the  open-loop  and  closed- 
loop  solutions  for  low-thrust  spacecraft  maneuvers.  As  shown  in  Tables  1,2 ,t;  increases  by  exerting  disturbances  on  the 
system  and  it  can  be  computed  in  this  method.  Also,  variations  of  control  value  at  the  end  of  the  trajectory  are  determined. 
Finally,  the  results  illustrate  the  viability  of  the  NOC  method  as  analytical  closed-loop  guidance  laws  for  low-thrust  orbital 
maneuvers  regarding  the  robustness  of  these  guidance  methods  against  disturbances. 
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